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Pseudo-Newtonian gravitational potential describing the gravitational field of static and 
spherically symmetric black holes in the universe with a repulsive cosmological constant 
is introduced. In order to demonstrate the accuracy of the pseudo-Newtonian approach, 
the related effective potential for test-particle motion is constructed and compared with 
its general relativistic counterpart given by the Schwarzschild-de Sitter geometry. The 
results indicate that such an approach could be useful in applications of developed New- 
tonian theories of accretion discs in astrophysically interesting situations in large galactic 
structures for the Schwarzschild-de Sitter spacetimes with the cosmological parameter 
y = (1/3)AM 2 < 10" 6 . 

Keywords: Pseudo-Newtonian approach; relativistic approach; gravitational potential; 
effective potential; circular orbit; energy; cosmological parameter; black hole 

1. Introduction 

Recent cosmological observations of distant supernova la type expl osio ns suggest 
an accelerating universe; starting at the cosmological redshift z ~ the accel- 
erated expansion should be generated by the so-called dark energy. These results 
are in accord with detailed observations of the anisotropies of the microwave cos- 
mic background radiation indicating that the dark energy represents ~ 73% of the 
energy content in the observable univers the sum of energy densities is very close 
to the critical energy density p cr it correspondin g th e parabolic universe in concor- 
dance with the inflationary cosmological modeP^. There is a variety of possible 
candidates for the dark energjE^. Most of the related models have the equation of 
state with parameter w = p/ p varying during the cosmological expansion. 

Nevertheless, the recent observational data indicate that the allowed equation 
of state is very close to the case of constant w = — 1 corresponding to the repulsive 
cosmological constant A > with pjy ~ 0.73p C rit • Therefore, it is quite important to 
consider the influence of A > in astrophysically relevant situations quite seriously. 



We have tested possible effects of A > in astrophysically motivated problems-*-', 
investigating namely the properties of black-hole solutions of Einstein equations for 
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the test-particle motion in Schwarzschild-de Sitter (SdSjSland Kerr-de Sitter (KdS) 9 
spacetimes and the test perfect fluid in spacetimes. Further, the 

test-particle and fl uid pr operties were treated also in the framework of the optical 
reference geometry ^^ ^ allow ing; introduction of inertial forces in the intuitively 
natural Newtonian wajff^HSl 

The influence of the repulsive cosmological constant on the black-hole spacetimc 
structure can be properly represented by the dimensionless cosmological parameter 
y = (A/3)M 2 . For SdS black holes admitting existence of stable circular geodesies, 
i.e., existence of accretion discs, the cosmological parameter y < y ms , e = 0.000237 
El. The cosmological tests using the supernova magnitude-redshift r elati on and the 
Cosmic Microwave Background Radiation fluctuations measurements^^ imply A w 
10~ 56 cm~ 2 , and thus very low values of y for astrophysically realistic black holes. 
In fact, y ~ 10~ 24 for super-massive black holeJISl strong optical observable effects 
could be expected for super-giant black holea^ (or clusters of galaxies) with M > 
10 15 M Q 0. 

For astrophysically realistic SdS black holes, the strong gravity near the black 
hole horizon ~ 2M weakens with distance growing and at r >> M can be de- 
scribed quite well by the Newtonian theory. However, the Newtonian theory looses 
its validity near the so-called static radius r s ~ y~ x / 5 M, where the repulsive effect 
of the cosmological constant starts to be relevant up to the other strong gravita- 
tion region near the cosmological horizon r c ~ y~ x l 2 M. Therefore, the cosmological 
constant has relevant influence on the structure of disc configurations (both geo- 
metrically thin and thick) introducing quite naturally outer edge of the accretion 
discs^lM 

It should be stressed that all of the relevant effects of A > on the black-hole 
structure are quite well expressed in the SdS spacetimes, since the rotational effects 
of the black hole spin are concentrated into the region in close vicinity of the black- 
hole horizon, where the influence of A > can be abandoned for realistic values of 
the black hole mass and the relict cosmological constant. Of course, the efficiency of 
the accretion processes is controlled by the rotational effects in the innermost parts 
of the disc, where the KdS spacetime structure is relevant!^. For near-extreme black 
holes the efficiency takes large values ~ 0.4 as compared with the SdS spacetimes 
where the efficiency ~ 0.059. Therefore, in studying the large scale properties of disc 
structures, investigation of the SdS spacetime is quite sufficient, only the accretion 
efficiency has to be given by the K(dS) spacetime structure, governing the innermost 
parts of the disc. 

For this reason, it is worth to realize more detailed studies of the disc structures 



a For primordial black holes in the very early universe, with expected high values of effective 
cosmological constant, the values of y can be much closer to 2/ m s,e< Considering the electroweak 
phase transition at T ew ~ lOOGeV, we obtain an estimate of the primordial effective cosmological 
constant A wc ~ 0.028cm~ 2 , while at the level of the quark c onfine ment at T qc ~ lGeV we obtain 
A qc ~ 2.8 X 10 — 10 cm - 2 and consequently higher values of yEES 
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around supermassive black holes reflecting the influence of the repulsive cosmologi- 
cal constant. There is one especially important problem that could be hardly solved 
in the framework of full general relativistic approach, namely influence of A > on 
the structure of self-graviting disc. We expect that the pseudo-Newtonian approach 
could be successful. For such purposes, we introduce here a pseudo-Newtonian 
potential tp of the SdS spacetimes that could enable us to use directly standard 
techniques developed in the framework of Newtonian physics, e.g., the Newtonian 
discoseismology^, taking the influence of the repulsive cosmological constant and 
some relativistic effects into account. 

There is a variety of different approaches in defining the pseudo-Newtonian 
gravitati onal potential describin g black holes and various aspects of their spacetime 

II 811 QI90I91 1 00 9^ OA 9 EH ioi I 

structur c! 10 ' l: ^ u ' zl ' ^ ' ^ J '. In the case of Schwarzschild spacetimes, it seema^ 
that to reflect the accretion disc properties, the most convenient is the Paczynski- 
Wiita (P-W) gravitational potentiate tpP-w = - 2GM/c 2 ). It enables us to 

calculate positions of the marginally stable and bound circular orbits at the same 
radii as follow from the general relativistic calculations. 

Originally, the P-W potential was introduced by a guess, when attempting to 
include the Schwarzschild radius into the Newtonian gravity. However, there is a sim- 
ple heuristic method for derivation of the pseudo-Newtonian potentials that yields 
the P-W potential. The same method was used for the derivation of the pseudo- 
Newtonian gravitational potential for the equatorial plane of Kerr spacetimes as 
welP^l. Then the position of the marginally stable circular orbit corresponds to the 
position determined by using the general relativistic approach, and differences in 
positions of marginally bound circular orbit determined in both the ways are rela- 
tively small. In this paper, this heuristic method is just used for derivation of the 
pseudo-Newtonian gravitational potential for the SdS spacetimes. 

Defining such potential, we have to reflect properly both the gravitational at- 
traction of the black hole and the repulsive effects of the cosmological constant. The 
P-W potential, describing with high precision properties of accretion processes in 
the field of Schwarzschild black holes has naturally defined behaviour at infinity. In 
the SdS spacetimes, we have to find a proper reference point that could in a well 
define d way serve in a similar sense as infinity in asymptotically flat spacetimes. It is 



the gravitational attraction of black hole acting on matter is exactly balanced by the 
cosmic repulsion, i.e., test-particles feel no force, similarly to asymptotic infinity in 
flat spacetimes. Concentrating on properties of stationary disc configurations in the 
SdS spacetimes, the static radius is appropriately chosen; we do not consider here 
properties of the spacetime near the cosmological horizon of the SdS spacetimes, 
nevertheless, the pseudo-Newtonian potential we are using here reflects the basic 
spacetime property near the cosmological horizon, since it diverges there. 

In the pseudo-Newtonian potential we properly fix position of the marginally 
stable orbit, marginally bound orbit, the static radius, and we establish the particle 
energy at these radii in agreement with the exact general relativistic values. 




that such a role could be attributed to the so-called static radius, whe: 



■re 
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2. SdS spacetime and its pseudo-Newtonian gravitational potential 

In the standard Schwarzschild coordinates (t,r,9,cj>), and the geometric system of 
units (c = G = 1), the SdS spacetimes are determined by the line element!^ 

2M A 9 \ 9 
= [ 1- — - -r 2 J dt 2 

l_^_| r 2^ 1 dr 2 + r 2 ^2 + ^2 M ^ ) (1) 

where M is the mass parameter of the spacetimes. It is useful to introduce the 
dimensionless parameter y = AM 2 /3 and use dimensionless coordinates t — > t/M, 
r — > r/M, which is equivalent to putting M = 1. 

Singularities of the line element, i.e., black-hole and cosmological horizons, are 
given by the relation 1 — 2/r — yr 2 = 0, thus by solutions of the equation 

r — 2 

y = y* = —> ( 2 ) 

which can be expressed in the form 

2 7T + £ 2 7T-£ 

rv, = ; cos , r c = — == cos , ( d ) 

where ^ = cos -1 (3^/3^). Both the horizons exist for < y < y c rit = 1/27, separating 
the spacetimes into two dynamic and one static regions. For y = y cr ;t, both the 
horizons coalesce at the radius of the photon circular geodesic at r = r p h = 3. For 
V > 2/criti the horizons disappear and the SdS spacetimes become dynamic naked- 
singularity spacetimes (see Fig.Q]). 

The heuristic method^ enabling us to define the P-W gravitational potential 
and the pseudo-Newtonian gravitational potential for the equatorial plane of the 
Kerr spacetimes is based on the knowledge of exact general relativistic relations for 
the angular momentum per particle mass L c and energy per particle mass E c of 
the circular geodesies. The Newtonian gravitational potential V'n for central grav- 
itational fields is related to the Newtonian angular momentum per particle mass 
Zn.c of the circular orbits by the relation dip^/dr = l 2 ^ c /r 3 . The pseudo-Newtonian 
gravitational potential ip is defined by the Newtonian relation with the transposition 
Zn.c — ► L c /E c = Z c 0, i.e., by the relation 

Note that the described method of pseudo-Newtonian determination works quite 
well in spherically symmetric (non-rotating) spacetimes, or in the equatorial plane 
of axially symmetric (rotating, e.g. Kerr or KdS) spacetimes. However, it is much 



b The quantity l c = L c /E c plays its role only when the pseudo-Newtonian (e.g, P-W) gravitational 
potential is denned. Later, standard Newtonian quantities in Newtonian theory are used along 
with the pseudo-Newtonian gravitational potential. 
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Fig. 1. Functions ?/h(r) (solid) determining radii of horizons, y m h (dashed) determining 
marginally bound circular orbits, and y,ms(r) (dotted) determining marginally stable circular or- 
bits. Regions above the solid curve correspond to dynamic parts of SdS spacetimes. 



more complicated task to find a pseudo-Newtonian potential for the regions outside 
the equatorial plane of the rotating spacetimes, because of a non-trivial influence 
of the dragging of inertial frames. There is a need to upgrade this method^ or use 
completely different way of the gravitational potential definitional 

In general relativity, the circular orbits in the SdS spacetimes correspond to 
extrema of the effective potential 



V e s(r;L,y) = 



1 yr' 

r 



1 + — 



1/2 



-1/2 



-1/2 



(5) 

(6) 
(7) 

2(2 -r + r 3 y) ' ^ 
with an integrating constant K. We calibrate the potential ([8]) by using the condition 



whereas the related constant of motion take the form 
L c (r;y)^[r(l-yr 3 )} 1/2 (l - * 

E c {r nj )^(l-- r -yA fl- 3 - 



Thus, the gravitational potential ((4]) can be written in the form 



ip(r = r s ) = 0, r s 



-1/3 



(9) 



where r s is the static radius defined above; the gravitational f orce, defined in the 
general relativistic framework of the optical reference geometry^^Sl disappears at 
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Fig. 2. Pseudo-Newtonian gravitational potential ip(r; y) for three values of the cosmological 
parameter y = (dotted), y = 10~ 5 (dashed) and y = 10~ 4 (solid). The values of y are chosen 
high enough to exhibit in most illustrative way the structure of SdS spacetimes allowing the 
existence of accretion discs. 



r = r s . Existence of equilibrium position of free particles with covariant energy 
representing the rest energy 

E s (y) = E c (r = r s :y)= (l-'Sy^y 12 (10) 

is possible there. Therefore, it is natural to consider the static radius as the 'starting 
point', corresponding to infinity in asymptotically flat spacetime. The proposed 
pseudo-Newtonian gravitational potential ([5]) then takes the form 

r 3 y ~ Sry 1 / 3 + 2 
^ (r; ^ 2(l-3,V3)( 2 - r + r 3, r (") 

diverges at the radii of horizons and vanishes at the static radius; for y = the 
potential tends asymptotically to zero (see Fig. [2]). 



3. Test-particle motion 

In the case of central gravitational fields, test-particle motion is confined to central 
planes and we can choose the equatorial plane for an individual particle. Following 
the Newtonian physics, the radial equation of the Keplerian equatorial motion can 
be written in the form 

2-UJ = e ~ V ^ (12) 

where e is the total pseudo-Newtonian energy per particle mass (energy hereafter) 
and v e g is the pseudo-Newtonian effective potential per particle mass (effective 
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potential in the following) defined by the standard relation 

I 2 

v e s(r;l,y) = ip + ^-j. (13) 

Here, ip(r; y) is the pseudo-Newtonian gravitational potential (fTTjl and I is the 
pseudo-Newtonian angular momentum per particle mass (angular momentum here- 
after) defined in Section O The circular Keplerian orbits (geodesies) correspond to 
the effective potential extremE^. Thus, their radii are governed by the function 



r 3 (l-r 3 y) 
(2-r + r 3 y) 2 



lt(r;y)= - y ■ »' (14) 



and the corresponding energies (effective potential extrema) are governed by the 
function 

6c(r; V) = 2(1-3^/3) - 2(2-r + r3y)2- (15) 
3.1. Marginally stable and marginally bound orbits 

Behaviour of the pseudo-Newtonian effective potential u e g qualitatively follows be- 
haviour of its general relativistic counterpart V e s (see Fig. We demonstrate this 
statement by comparing behaviour of the functions I 2 and e c , and their general 
relativistic counterparts L 2 and E c , defined by relations ©, 0, fT4")) and (fT5|) . 

The functions I 2 and e c diverge at the radius of the black-hole horizon and 
vanish at the static radius, while L 2 and E c diverge at the radius of the photon 
circular orbit, L 2 vanishes at the static radius, and E c does not (see Figs [4] and [5]) . 
Thus, the pseudo-Newtonian approach implies that existence of circular orbits is 
limited by the static radius from above, in accord with the relativistic approach, and 
by the black-hole horizon from below, in contrast to the relativistic approach, where 
the lower limit is given by the radius of photon orbit. This is caused by the fact that 
we do not obtain the photon circular orbit in the pseudo-Newtonian approacrj^. 

As well as in the general relativistic approach, stable circular orbits, determined 
by the minima of the effective potential v g, satisfy the condition d r l 2 > 0. The 
unstable circular orbits, determined by the maxima of v e s, satisfy the condition 



c Keplerian circular motion can be equivalently given also directly from the pseudo-Newtonian 
gravitational potential mil l using relations for orbital and angular velocitites, and for the angular 
momentum and energy 



However, in some sense, the method of effective potentiaJ23 combining the gravitational potential 
and potential of centrifugal forces, is more general, illustrative and convenient for our case. 
d The same situation occurs for the standard P-W pseudo-Newtonian potential denned for the 
Schwarzschild spacetimeii&l 
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Fig. 3. Pseudo-Newtonian v e ff(r; I, y) (thick) and relativistic V e g (r; L, y) — E s (y) (thin) effective 
potentials for test-particle motion with the pseudo-Newtonian I 2 = l 2 (r = 7) = 13.4) and rela- 
tivistic L 2 = L 2 (t = 7) = 11.8 angular momenta in the SdS spacetime with y = 10 — 4 . Both the 
potentials determine stable circular orbits at r = 7. The dotted lines suggest determination of 
the energy differencies (potential barriers) between the outer unstable orbits, determined in the 
pseudo-Newtonian and relativistic ways, and the fixed central stable orbit at r = 7. 



d r lc < 0. Inner and outer marginally stable circular orbits correspond to extrema 



Z„ lsi (y) and i^iso(y) °f tne function l\. They are governed by the function 



Vms{r) 



6 



r 3 (4r - 15) ' 



(16) 



just as in the general relativistic approach, because loci of extrema of the functions 
ll and L\ coincide. The stable circular orbits are then limited by the condition 
4yr 4 - 15yr 3 - r + 6 < 0. 

The function y ms vanishes at the radius r = 6, corresponding to the marginally 
stable circular orbit in the Schwarzschild spacetime, diverges at r = and at r — 
15/4, and has minimum at r = r p h, where y m s(f = fph) = J/crit- The function is 
irrelevant at the range < r < 15/4, where y ms > j/ C rit- The physically relevant 
part of y ms is located at r > 6. Its maximum is located at r ms o = 15/2 and the 
corresponding maximum is y ms (f = ?"ms,e) = J/ms,e = 12/15 4 (see Fig.[l]). This value 
represents the limiting value for the SdS spacetimes admitting stable geodesies, i.e., 



accretion discs^ 1 . Behaviour of the effective potential determining the innermost 
and outermost stable circular orbits is illustrated in Fig. |6] 

The marginally bound circular orbits (inner and outer), corresponding to two 
unstable circular orbits with the same energy and appropriately chosen angular 
momentum, are determined by the condition 

e c (r = r mhi ; y) = e c (r = r mbo ; v) = e mh (y). (17) 

The related radii r m bi and r m b Q can be determined by using the following numerical 
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Fig. 4. Radial profiles of pseudo-Newtonian l^(r;y) (thick) and relativistic L^(r;y) angular mo- 
menta determining circular orbits (effective potential extrema) in SdS spacetimes with y = 10 — 4 . 
We show an example of typical behaviour of the angular momenta, each with two extrema. The 
solid vertical line denotes the black-hole horizon position; the dashed vertical lines denote pho- 
ton and marginally stable circular orbits positions; the dashed-dotted vertical line denotes the 
static radius position. The dotted lines suggest determination of the difference in positions of the 
outer unstable circular orbits, determined in the pseudo-Newtonian and relativistic ways for fixed 
position of the central stable orbit at r = 7. 

procedure. From equation (|14p . we express two radii ri(l 2 ;y) and rn(l 2 ;y) corre- 
sponding to the inner and outer unstable circular orbits (marginally bound orbits 
radii candidates). Taking, moreover, the condition ()17ll into account, i.e., solving 
the equation e c (r = rj;y) — e c (r = rjj;y), where e c (r;y) is given by relation (|15[) . 
and considering formula (fT7|) once more, we obtain function y m b(V) implicitly de- 
termining the radii of the marginally stable orbits r m u(y) and r mrjo (j/) (see Fig. [7|)- 
General relativistic formulas imply the same results, because 

e c {r-y) = \[E^{y) - E- 2 (r-y)] , (18) 

and for fixed values of r, corresponding to r m bi and r m bo, there is E c (r = r m bi) = 
E c (r = r mbo ), and thus e c (r = r mb i) = e c (r = r mbo )- 

3.2. Potential barriers 

For values of the cosmological parameter y allowing existence of the stable circular 
orbits, i.e., < y < y ms , C j the functions l^(r;y) and L^(r;y) have two extrema 
governed by the function y ms (r) (see Fig. |4]). These extrema correspond to the in- 
nermost and outermost stable circular geodesies with I 2 = I 2 - and I 2 = l 2 a „, and 

1 Illbl IIloU i 

L 2 = i, 2 nsi and L 2 — £, 2 nso , respectively, which determine extension and accretion 
efficiency of thin, Keplerian discs. The marginally stable perfect fluid configurations 
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Fig. 5. Pseudo-Newtonian e c (r; y) (thick) and relativistic E c (r; y) — E s (y) energies of circular or- 
bits in SdS spacetimes with y = 10 — 4 . The dotted lines suggest determination of the energy differ- 
ences (potential barriers) between the outer unstable orbits, determined in the pseudo-Newtonian 
and relativistic ways, and the central stable orbit at r = 7 (for the dashed-dotted and dashed 
vertical lines meaning see Fig. [2). 

with l 2 (r, 9) =const and I 2 S (^ s ;, ^ so ) are determined by intersections oil 2 =const 
and I 2 curve that determine the center and edges of the fluid configuration^^. There- 
fore, we realize a detailed comparison of the potentials u c ff and V e g, concentrating 
on the differences in positions of the inner, central and outer circular orbits (deter- 
mined in the pseudo-Newtonian and relativistic ways), and the differences in the 
corresponding values of both the potential barriers. 

We start with analyzing the pseudo-Newtonian and relativistic energy differences 
between the outer and inner marginally stable circular orbits, which determine 
efficiency of the accretion processes in Keplerian discs 

Ae ms (?/) — ^mso ^msi; (-^) 

AE ms (y) = E mso - E msi . (20) 

The energy differences Ae^s an d AE ms are compared in Fig. [8] Clearly, for astro- 
physically relevant value d I I of y < 10~ 15 , the pseudo-Newtonian and relativistic 
efficiencies are very close. The differences of the pseudo-Newtonian and relativistic 
efficiency, can be characterized by the quantity 

x(y) = Ae ™ ~ AE ™ x ioo%, (2i) 

the dependence of which on the cosmological parameter is illustrated in Fig. [H] In 
the limit of y — > 0, i.e., in the Schwarzschild spacetimes, there is Ae ms = 0.063, 
AE ms = 0.057 and x = 9.3%, and these values hold for the astrophysically relevant 
values of y. 
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Fig. 6. Marginally stable circular orbits. Pseudo-Newtonian v e g(r; I, y) (thick) and relativistic 
V e fj(r; L, y) — E B (y) effective potentials for test-particle motion with the pseudo-Newtonian I 2 = 
l 2 (r = 6.2) = 13.3 (upper figure), I 2 = l 2 (r = 12.3) = 14.8 (lower figure), and relativistic L 2 = 
L 2 (r = 6.2) = 11.7 (upper figure), L 2 = L 2 (r = 12.3) = 13.2 (lower figure) angular momenta in 
the SdS spacetime with y = 10 -4 . Positions of the marginally stable circular orbits correspond to 
the positions of the potential inficx points. 

Further, it is of astrophysical relevance to compare dependence of the pseudo- 
Newtonian and relativistic characteristics of marginally stable discs with uniform 
distribution of angular momentum l 2 (r, 9) =const. Assuming a fixed radius of stable 
circular orbit corresponding to the disc centre, we can determine and compare the 
outer edge of the disc and the potential energy well for the matter of the disc. 

By fixing a value of the radius of stable circular orbit r = r st , and calculat- 
ing the corresponding values of both the pseudo-Newtonian and relativistic angular 
momenta l\ (r = r st ; y) and L\ (r = r st ; y) (see Fig. |4]) , we can determine the radii 
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Fig. 7. Marginally bound circular orbits. Pseudo-Newtonian v e g (r; I, y) (thick) and relativistic 
V e g(r; L,y) — E s {y) effective potentials for test-particle motion with the pseudo-Newtonian I 2 = 
l 2 (r = 15.3) = 13.7 and relativistic L 2 = L 2 (r = 15.3) = 12.2 angular momenta in the SdS 
spacetime with y = 10~ 4 . Positions of the marginally bound circular orbits correspond to positions 
of effective-potential maxima. 




Fig. 8. Left: Efficiency of the Keplerian discs accretion determined in the pseudo-Newtonian 
Ae ms ({/) (thick) and relativistic AEms(y) ways. Right: Characteristic quantity \ describing the 
inaccuracy in the determination of the Keplerian accretion efficiency. 



of the corresponding outer unstable circular geodesies (edge of the disc) r uno (r st ; y) 
(calculated by using 1%) and R U no{r s t',y) (calculated by using L%), and their differ- 
ence 



Ar(r st ;y) = r uno - R uno . 



(22) 
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The inaccuracies in the pseudo-Newtonian determination of the position of the outer 
unstable orbit, in dependence on the position of the stable central circular orbit, 
can be characterized by the quantity 

V) = p ^ x 100%. (23) 

We also calculate the pseudo-Newtonian energies at the radius of the stable circular 
orbit e s t(r s t; y) = e c (r = r st ; y) and at the related outer unstable orbit e uno (r st ; y) = 
e c (r = r uno ;y), and their difference (potential barrier) 

Ae(r st ; y) = e uno - e st . (24) 

The same can be done in the relativistic case, i.e., we calculate the energies 
E s t(r st ;y) = E c (r = r st ;y) and E uno (r st ;y) = E c (r = R uno ;y), and their differ- 
ence (potential barrier) 

AE(r st ;y)=E uno -E st . (25) 

The inaccuracies in the pseudo-Newtonian determination of such a potential barrier, 
in dependence on the position of the stable central circular orbit, can be character- 
ized by the quantity 

Dependencies of both the characteristic quantities £ and rj are illustrated in Fig. [5] 
for selected values of y. 

From the astrophysical point of view, it is better to study the differences be- 
tween the pseudo-Newtonian and relativistic determination of the radii of the outer 
unstable circular orbits and potential barriers between the inner and outer unstable 
orbits, when the position of the inner unstable circular orbit at r — r un ; is fixed. In 
analogy with the previous case, we define the characteristic quantities 

arum; y) = -5-^ x 100%, (27) 

-'•'uno ' uni 

fj(rnm;y)= Ae A ^ x 100%, (28) 

where 

Af(r uni ; y) = f uno - R uno , (29) 

Ae(r uni ; y) = e uno - e uni , (30) 

AE(r uni ;y) = E uno - E uni . (31) 

The radii f U no(?"uni; y) and R U no(r un i; y) are now determined by calculating l%{r = 
r uni ; y) and L\(r = r u „i;?/), respectively. Corresponding energies are e uno (r uni ;y) = 
e c (r = r uno ;y) and E uno (r uni ;y) = E c (r = i? uno ;y), and there is e U ni(r uni ;y) = 
e c (r = r unl ;y) and E un i(r unl ;y) = E c (r — r unl ;y). Dependencies of both the char- 
acteristic quantities ^ and f\ are given in Fig. [10] for selected values of y. 
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Fig. 9. Left: Characteristic quantity £(r s t;y) describing the inaccuracies in the pseudo- 
Newtonian determination of the positions of the outer unstable orbit; Right: Characteristic quan- 
tity n(r B t;y) describing the inaccuracies in the pseudo-Newtonian determination of the potential 
barrier between the outer unstable and central stable circular orbits in the SdS spacetimes with 
y = 1CT 4 (solid), y = 10~ 5 (dashed-dotted) y = 10~ 6 (dashed) y = 10~ 7 (dotted) y = (thick 
dotted) when position of the stable circular orbit is fixed. The values of y are chosen high enough 
to exhibit in most illustrative way the influence of y on the calculated quantities. 




34567 34567 



Fig. 10. Left: Characteristic quantity £(r unl ;y) describing the inaccuracies in the pseudo- 
Newtonian determination of the positions of the outer unstable orbit; Right: Characteristic quan- 
tity v( r nni',y) describing the inaccuracies in the pseudo-Newtonian determination of the potential 
barrier between the outer unstable and inner unstable circular orbits in the SdS spacetimes with 
y = 10~ 4 (solid), y = 10~ 5 (dashed-dotted) y = lO" 6 (dashed) y = 10~ 7 (dotted) y = (thick 
dotted), when position of the inner unstable circular orbit is fixed. 

Clearly, the precision in the determination of the radii of the unstable circular 
orbits and potential barriers, characterized by the quantities £ (£) and r\ (fj) is high 
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enough, when y < 10~ 6 . Note that for y — ► 0, the positions of unstable circular 
orbits r uno (f uno ) — > oo and i? U no(^uno) - * 00 , and thus — > 0. In fact, it is 
irrelevant to study £ (£) for y = 0, because there arc no outer unstable orbits 
in the Schwarzschild spacetime. On the other hand, there is e uno (e uno ) — ► and 
E uno (E uno ) — » 1 for y — » and r uno (f\mo) — » oo and R un o(Runo) — » oo. Thus, 
in the case 3/ = 0, 77 (77) describes the differences in the pseudo-Newtonian and 
relativistic determinations of the potential barrier between the stable (unstable) 
orbit and infinity. 

4. Conclusions 

We have shown that the gravitational field of spherically symmetric and static 
black holes in the universe with a positive cosmological constant, described by 
the Schwarzschild-de Sitter solution of Einstein's equations, can be alternatively 
described with a relatively high precision by using appropriately defined pseudo- 
Newtonian gravitational potential. 

The presented gravitational potential satisfies important conditions. It admits 
existence of the static radius, diverges at horizons, and it gives marginally stable 
and bound orbits at radii exactly equal to those given in the relativistic expressions. 
The energy difference of these orbits is close to the relativistic relations. 

We have tested the potential correctness by comparing some pseudo-Newtonian 
results concerning the test-particle gcodctical motion with the general relativis- 
tic ones. We have been interested in astrophysically relevant situations. Thus, we 
have chosen range of the cosmological parameter (0 < y < 0.000237), allowing the 
existence of stable circular geodesies, i.e., existence of accretion discs. The differ- 
ences in the pseudo-Newtonian and relativistic calculations of the accretion pro- 
cesses efficiency have been also studied, whereas the related inaccuracies of the 
pseudo-Newtonian efficiency determination are less then 12% for y < I0~ 6 , and 
in astrophysically realistic cases of y < 10 -15 , the inaccuracies are very close to 
those corresponding to the standard P-W potential, describing the Schwarzschild 
spacetime. 

Moreover, assuming the central circular orbit position fixed, we have compared 
results of the pseudo-Newtonian and relativistic calculations of the outer unstable 
circular orbits positions, and the related energy differences (potential barriers) be- 
tween the central stable (inner unstable) and outer unstable orbits. The comparison 
suggests that within a certain inaccuracy (less than 0.5% in the case of the radii 
determination, and less than 12% in the case of the potential barriers determina- 
tion, for y < I0~ 6 ), the pseudo-Newtonian results are in a good agreement with 
the relativistic ones. As for the potential barriers, we can state that the differences 
between the results of both the approaches tend to be smaller for y decreasing and 
approaching to y = 0. 

Note that, in principle the existence of three related different kinds of circular 
orbits, i.e., the inner unstable, central stable and outer unstable orbits, provides 
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us more alternatives for the pseudo-Newtonian gravitational potential testing pro- 
cedures than presented. We can fix the position of the outer unstable orbit and 
calculate the differencies in the determination of the central stable orbit and related 
potential barriers etc. We have not proceeded these comparisons, representing, in a 
way, alternatives of the presented ones. We prepare a detailed study of the appli- 
cation of the pseudo-Newtonian potential on the perfect fluid dynamics, extending 
thus further our intention to study accretion processes and the properties of ac- 
cretion discs in dependence on the cosmic parameter y in the pseudo-Newtonian 
way. 
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